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Il'ITRODUCTION 


The  subject  of  this  work  is  a study  of  the  ground  state  of  the 
Be  atom  by  the  configuration  interaction  method,  with  particular  em- 
phasis on  tha  total  energy  as  well  as  on  the  quality  of  the  wave  func- 
tion. The  techniques  used  here  can  be  applied  quite  in  general,  to  the 
calculation  of  energy  levels,  wave  functions  and  related  quantities  of 
isolated  atoms  in  their  bound  states. 

The  absolute  importance  of  finding  these  quantities  resides 
in  innumerable  applications  in  the  field  of  atomic  spectroscopy 
(energy  levels,  transition  probabilities,  hyperfine  structure),  togeth- 
er with  the  evaluation  of  average  values  giving  information  on  the 
shape  of  the  electronic  cloud,  and  of  all  the  physical  properties  of 
isolated  atoms  such  as,  for  example,  their  diamagnetic  susceptibility. 

The  relative  importance,  however,  which  is  primordial  in  any 
growing  field,  concerns  mainly  quantum  chemistry,  where  the  analysis 
of  accurate  v/ave  functions  offers:  a)  an  opportunity  for  testing  and 
improving  existing  theories  on  atomic  structure  and  perhaps  for  suggest- 
ing new  ones;  b)  a starting  point  for  better  atomic  calculations, 
and  c'  information  useful  in  molecular  calculations,  mainly  in  rela- 
tion to  the  search  of  appropriate  basic  functions. 

Ever  since  the  early  days  of  quantum  mechanics,  it  has  been 
realized  that  it  is  almost  impossible  to  solve  the  electronic  Schrd- 
dinger  equation,  and  people  have  concentrated  on  the  computational 
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point  of  viev;,  inventing  models  and  finding  methods  to  handle  the 
complicated  mathematics  involved  even  in  the  roughest  approximations. 

It  is  believed  that  some  intermediate  results  are  necessary  in  order 
to  make  possible  the  formulation  of  a many  electron  theory.  It  is  on 
these  intermediate  results  that  most  of  the  computational  work  is 
being  done  nowadays  and  ours  falls  in  this  category.  The  Hartree-Fock 
(KF)^-^^and  Hartree-Fock-Slater^^^  calculations  on  atomic  systems  are 
an  example,  and  they  are  those  which  have  been  best  studied  both  theo- 
retically and  from  the  practical  point  of  view. 

The  HF  approximation-})  or  independent  particle  model,  is 
regarded  by  many^^^  as  the  basis  for  the  understanding  of  the  electron- 
ic properties  of  atomic  and  molecular  systems  in  their  stationary 
states.  One  electron  properties  calculated  from  good  approximations 
to  the  HF  wave  functions  are  in  very  good  agreement  with  experiment 
(there  is  a theoretical  basis  to  account  for  this  fact^^^).  For  exam- 
ple, the  HF  electronic  densities  agree  very  well  with  X-ray  results^^^. 
Much  more  important,  the  basis  for  the  modern  picture  of  structural 
chemistry  in  terms  of  cores,  shells,  bonds  and  lone  pairs,  and  the 
treatment  of  excitation  processes  lies,  to  a great  extent,  in  the  or- 
bital description  of  these  phenomena.  To  v/hat  extent  these  pictures 
remain  correct  v/hen  going  to  approximations  more  accurate  than  HF 
is  one  of  the  interesting  problems  which  concerns  the  quantum  chemists 
of  today 

The  idea  that  molecules  are  made  of  atoms  which  retain  much  of 
their  individuality  after  bond  formation,  is  still  at  the  roots  of 
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any  molecular  theory,  and  it  has  been  verified  by  both  rigorous  and 
semiempirical  calculations  (the  hydrogen  atom  being  sometimes  an 
exception).  The  same  is  true  when  one  considers  the  removal  of  one 
electron  to  form  a closed  shell:  the  structure  of  the  closed  shell 
remains  very  much  the  same  after  ionization  of  the  outer  electron. 

Thus,  the  obtaining  of  very  accurate  atomic  wave  functions 
may  be  useful  to  investigate  bond  formation  and  the  simple  and  visual 
models  of  excitation  processes.  In  this  connection,  we  may  mention 
that  presently,  there  is  only  one  rigorous  method  to  calculate  binding 
energies  or  the  ionization  potential  of  an  outer  electron:  it  consists 
in  subtracting,  from  the  energy  of  the  final  state  (molecular  aggre- 
gate or  ionized  atom),  the  energy  of  the  corresponding  initial  state. 
Both  quantities  are  usually  so  large  compared  to  their  difference 
that  the  method  has  been  compared  to  weighing  the  captain  of  an  ocean 
liner  by  taking  the  difference  between  the  weight  of  the  ship  with 
and  without  captain.  Certainly,  the  derision  of  the  method  does  not 
eliminate  the  problem,  and  we  are  confident  that,  with  the  help  of 
the  new  generation  of  computers,  the  present  methods  of  calculation 
will  afford  very  accurate  wave  functions  and  binding  energies  of 
small  molecules.  This  will  help  us  to  realize  better,  what  the  magni- 
tudes of  the  relevant  quantities  involved  in  the  existing  theory 
(7) 

are  ' , as  well  as  to  improve  the  application  of  it  to  larger  systems. 

There  are  several  ways  of  calculating  atomic  structures,  of 
various  degrees  of  complexity  and  accuracy,  and  which  are  meant  to 
serve  different  purposes:  ours  is  to  get  as  close  as  possible  to  the 


exact  wave  function.  The  simplest  of  all  is  the  configuration  interac- 
tion (Cl)  method,  or  method  of  superposition  of  configurations.  Its 
simplicity,  together  with  some  pessimistic  reports  about  a convergence 
failure,  has  led  many  to  think  of  it  as  a "pedestrian  approach."  The 
truth  is  that  none  of  that  pessimism  is  well  founded  and  it  reveals 
mainly  the  deception  of  those  v/ho  have  attempted  to  make  use  of  this 
Tiethod  without  the  right  tools. 

In  this  work,  we  examine,  from  the  practical  point  of  view, 
some  of  the  v;ell-known  difficulties  of  the  method,  such  as  the  ques- 
tion of  how  to  select  the  basic  one  electron  functions  and  the  differ- 
ent configurations.  This  is  treated  in  Chapters  I and  II.  In  Chap- 
ter III  we  consider  our  calculation  step  by  step.  Finally,  in  Chap- 
ter 17  we  analize  our  results. 


CHAPTER  I 


SXPAKSIOli  OF  THE  WAVE  FUIICTION 

One  of  the  difficult  proDlems  in  quantum  mechanics  is  to 
solve  the  many  electron  Schrodinger  equation 

= E.  (1.1) 

1 11 

v;here  H is  the  many  electron  Hamiltonian  and  is  the  v;ave  function 

1 

of  an  electronic  state  with  energy  E^  . 

An  expansion  of  the  v;ave  function  is  the  starting  point  for 
any  theory  or  calculation  v/hich  concerns  Eq.  (1.1).  The  many  electron 
theory  of  atoms  and  molecules,  as  recently  developed  hy  Sinanoglu'^  , 
starts  from  the  so-called  cluster  expansion,  whose  advantages  are 
pointed  out  in  Pief.  (7):  it  is  a finite  expansion  for  a finite  num- 
ber of  electrons  and  it  possesses  such  a convenient  form  that  each 
term  can  he  investigated  and  approximated  in  the  framework  of  the  above 
mentioned  theory;  the  uniqueness  of  Sinanoglu*s  cluster  expansion 
is  determined  by  the  HF  orbitals  of  the  ground  state  configuration. 
Anuther  expansion,  the  configuration  interaction  expansion^^^ , con- 
sists of  an  infinite  sum  of  determinantal  functions  (Slater  determi- 
nonts)  made  up  of  a complete  set  of  one  electron  functions,  or  spin- 
orbitals.  In  the  Cl  method,  the  linear  coefficients  of  the  Cl  expansion 
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are  found  "by  the  Ritz  variational  procedure^^^ . Formlly,  this  method 
is  very  simple,  and,  also  formally,  it  can  he  easily  applied  to  many 
electron  systems;  on  the  other  hand,  there  are  some  unsolved  difficulties 
with  the  cluster  expansion  method  in  this  respect^^®^.  However,  hy 
the  method  of  succesive  partial  orthogonalization^'^^  any  v/ave  func- 
tion, in  particular  a Cl  v/ave  function,  can  he  transformed  uniquely 
into  the  cluster  expansion  form  and  analyzed  according  to  the  correspond- 
ing many  electron  theory. 

It  is  very  important  to  stress  that  all  these  expansions  are, 
in  principle,  merely  formal^ In  order  to  make  use  of  them,  it  is 
indispensable  that  they  converge  to  the  exact  solution  in  a way  which 
i practical  from  the  computational  point  of  view  and  which  is  physically 
meaningful.  (From  the  mathematical  point  of  view,  it  would  he  very 
useful  to  know  something  about  the  real  convergence  conditions,  hut 
unfortunately  the  problem  has  only  been  clarified  for  the  case  of  two 
electron  systems^^^^ . ) 


1.1.  Configuration  Interaction  Expansion 


The  basic  assumiption  is  the  existence  of  an  expansion  theorem 


Y(x) 


(1.2) 


for  every  normalizable  function  ^(x)  of  the  space-spin  coordinate 
x=(r,^)  of  a single  electron,  in  terms  of  an  orthonormal,  complete 
..nd  discrete  set  (fpx)j  of  one  electron  functions,  or  spin-orbitals. 
For  a many  electron  wave  function  ^ , Lowdin  applies  the  same  theo- 
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ren,  one  coordinate  at  a time^^^.  After  introducing  the  antisyrainetri- 
zer  (since  the  electrons  are  fermions)  he  obtains: 


i 


X,j ) 


K 


r'^1’X2»***’X>^t) 


C. 


K 


v/here 


(1.3) 


]5^^(x^  ,X^  > . . » »X^j)  — (i'll) 


Vi 


m 


is  a normalized  Slater  determinant  K corresponding  to  an  ordered  con- 
figuration defined  by  the  set  k<  l<m  .... 

The  expansion  (1.3) ,vathout  any  other  symmetry  restrictions 
than  to  be  antisymmetric  with  respect  to  exchange  of  any  two  space-spin 
coordinates,  allowed  LSwdin  to  treat  the  convergence  problem  from  a 
most  general  point  of  view^^^.  He  considered  the  first  order  density 
matrix^(x’]  x^) 


^(x^  j x^ ) — (x^  jXgjX^, . . . ,Xj^)  ^ (x^  jXgjX^, . . . ,x,j)  dXgdx^. . .dXjj 


k,l 


and  observed  that  the  diagonal  elements  were  given  by 

1^/ 


whore  K runs  over  all  configurations  containing  the  specific  index  k. 
Then,  a new  basis  » such  that  ,was  introduced, 
so  that  ^ (x*  I x^)  became  diagonal: 


I (X’  I X^)  =ZX(x*)X  (x^)  nj^ 


(1.4) 
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The  's  associated  with  the  exact  wave  function  were  called  natural 
s'pin-orhitals  (the  others  are  called  approximate  natural  orhitals; 
both  have  symmetry  properties  if  these  are  included  in  ^(13). 

Unless 

explicitly  stated,  we  use  the  words  natural  spin-orhitals  (NSO)  for 
both  kinds  of  functions).  Vfnen  the  expansion  (1.3)  was  transformed 
into  the  j basis,  a long  forgotten  separation  theorem  due  to  Cau- 
chy^^^^  could  be  applied,  and  the  new  expansion  coefficients  were 
found  to  satisfy 


r 

k=l 


[kl  2 
lA 


Kl  >/ 


t f: 

'iC=l  K 


'K 


(1.5) 


the  index  k on  the  left  running  over  the  r natural  orbitals  associat- 
ed with  the  highest  occupation  number  n^  (see  Eq.  (1.4)),  and  that 
on  the  right,  running  over  any  r spin-orbitals  of  This  equation 

is  valid  for  the  exact  v/ave  function  as  well  as  for  any  approximation 
to  it.  Lowdin  concludes  by  pointing  out  that  Eq.  (1.5)  expresses  an 
optimum  convergence  property  of  the  natural  expansion. 


1.2.  Practical  Imuortance  of  ESO  Exuansions 


Here  we  are  going  to  consider  Eq.  (1.5)  more  in  detail,  so 
as  to  specify  under  v;hich  additional  conditions  it  expresses  an 
optimum  convergence  property.  Let  us  consider  the  natural  expansion 

and  let  us  subdivide  the  Slater  determinants  (SD)  into  two  classes:  L 
anc.  L’.(A  justification  for  this  will  be  found  at  the  end  of  Section 
(3.2).)  Those  in  class  L are  composed  of  natural  orbitals  with 
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occupancy  nuater  n^>  (see  Eq.  (1.4));  r>H.  Those  in 

class  L’  my  have  some  SD's  (an  infinite  number  perhaps)  which 
include  those  IISO's  showing  up  in  SD's  of  class  L,  hut  in  addition, 
they  may  have  other  ESO's  with  n^  4 • Equation  (1.5)  roay  now 

he  expressed  as 


or 


N J 


all  L 


■±  ' 

k=l  L'  ^ 


,|  > 'flu 

Cl 

‘ all  L 

2 _J1 

+ 


k=l  L 


(1.6) 


If  a very  careful  selection  of  L- type  SD's  is  made,  the  second  terms 
on  both  the  left  and  right  hand  sides  of  Eq.  (1.6)  turn  out  to  he 
extremely  small  and 

2 


I 


all  L 


1a 

^ > 

c. 

* L 

~ all  L 

L 

(1.7) 


expresses  the  optimum  convergence  property  in  a sense  v/hich  is  mathe- 
matically understandable.  For  two  electron  systems,  the  second  term 
on  the  left  vanishes  and  Eq.  (I.7)  is  mandatory  under  any  circum- 
stances. For  larger  systems,  however,  the  necessary  conditions  for  the 
fulfillment  of  Eq.  (1.7)  are  unknown.  Not  knowing  other  kinds  of 
linear  transformations  of  the  basic  spin-orbitals  leading  to  equations 
similar  to  (1.4)  or  (1.6),  i.e.,  leading  to  equations  v/hich  under 
som^  favorable  circumstances  may  conform  to  Eg.  (1.7) « it  is  worth- 
while to  examine  the  consequences  of  using  NSO's  in  practical  appli- 
cai;lons.  It  is  with  this  philosophy  that  we  undertake  to  experiment 


with  NSO  expansions. 
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We  Relieve  it  is  important  that  the  starting  wave  function  he, 
in  a hroad  sense,  a very  good  approximation  to  the  exact  solution.  It 
may  he  thought  that,  as  the  exact  solution  is  approached,  the  approxi- 
mate C.I.  v/ave  function  becomes  less  dependent  on  a linear  transfor- 
mation of  its  spin-orhital  basis,  and  that  thus,  Eq.  (1.7)  tends  to- 
wards an  equality.  This  is  incorrect.  The  answer  obviously  depends 
on  the  relationship  betv;een  the  starting  and  the  final  one  electron 
basis,  and  this  cannot  be  knov/n  unless  the  calculation  is  performed. 

1.3.  Application  of  NSO  ideas  to  Be  Ground  State 

Let  us  return  now  to  the  tv/o  electron  problem:  Lbwdin  and 
Shull  considered  the  mean  quadratic  deviation  between  the  exact 
solution  and  a truncated  C.I.  expansion  of  a two  electron  wave  func- 
tion^^^^  , and  proved  that  this  quantity  is  a minimum  for  the  case  of 
the  natural  expansion,  when  the  latter  is  truncated  to  any  order  r^^^^  . 

/T 

As  these  authors  pointed  out,  the  preceding  situation  is  equivalent 
to  the  diagonalization  of  a quadratic  form.  We  shall  use  this  fact 
in  the  following  context:  a four  electron  separated-electron-pair 
v;ave  * unction  n 

f\  =/4  [k(1,2)  L(3,4)] 

has  been  shov;n  to  be  a very  good  approximation  for  the  ground  state 
of  . ’We  assume  then,  that  an  essentially  correct  K(l,2)  func- 

tion which  is  presumably  localized  in  the  K shell,  could  be  obtained 
by  letting  the  L shell  v/ave  function  be  (2s)^  , i.e.  a doubly  occu- 
p*.iU  2s  orbital.  We  expand  K(l,2)  as 
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K(l,2)  =/4^_  ^ 


k<l 


kl 


(1.8) 


are 


in  vMch,  of  course,  has  opposite  spin  to  Lj/  ^ and  they 

always  conpled  to  give  a S two  electron  function.  The  natural  K 

shell  spin-orhitals  are  obtained  by  diagonalizing  the  quadratic  form 
in  Eq.  (1.8),  and  since  they  correspond  to  a two  electron  function, 
the  truncated  natural  expansion  of  K(l,2)  will  have  optimum  properties 
minimum  mean  quadratic  deviation  with  respect  to  the  original  K(l,2), 
or,  what  is  equivalent,  maximum  overlap  with  it.  In  this  way,  a large 
basis  can  be  effectively  reduced,  and  what  is  important  here  is  not 
the  obtaining  of  a basis  which  represents  K(l,2)  of  Eq.  (1.8)  in  an 
optimum  way,  but  rather,  to  get  a good  basis  in  which  to  expand 
that  localized  portion  of  the  total  wave  function,  the  "physical" 

K sh-.ll,  in  the  general  framework  of  an  unrestricted  Cl  calculation. 

Since  all  the  basis  orbitals  are  orthogonal,  the  natural  spin- 
orbitals  of  K(l,2)  are  identical  with  those  of  0.^^^ 


_Q(K)  ^K(l,2)  (2s)^} 


(1.9) 


except  for  the  addition  of  two  2s  spin  orbitals  with  occupancy  number 
equal  to  one. 

'«e  perform  the  same  treatment  with  a wave  function  n 


(1.10) 


in  tnis  fashion,  we  obtain  much  information  about  the  best  selec- 
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tion  of  "basis  spin-or"bitals  for  our  Cl  calculation.  The  most  impor- 
tant one  is  to  get  an  answer  to  the  question  of  how  many  spin-or"bitals 
of  a given  symmetry  type  are  energetically  effective  in  descri"bing 
a given  shell.  V/e  reserve  further  comments  for  Chapter  III. 


CIIAPTER  II 


SIjlKENTS  OF  A CONFIGURATION  INTERACTION  TREATI-IENT 


2.1.  The  Hamiltonian 

For  small  atoms,  the  relativistic  effects  are  small^^®^  and 
is  generally  accepted  that  all  aspects  of  atomic  "behaviour,  with 
the  exception  of  hyperfine  structure,  can  still  "be  accounted  for  if 
one  neglects  them.  In  doing  so,  the  hamiltonian  is  greatly  simplified. 

The  non-relativistic  Hamiltonian  H^,  in  a frame  of  reference 
(x^,X2>x^)  centered  at  the  nucleus,  is  given  "by 


^T  ^fn 


jn 

2K 


Vx,  *Vx 


(2.1) 


in  atomic  units  (a.u.;  1 a.u.  = twice  the  non-relativistic  energy  of 
the  "infinite  mass"  hydrogen  atom  ground  state).  is  the  conven- 
tional "fixed  nucleus"  non-relativistic  Hamiltonian 


-i).>  i 

^ 2 ^ r.  -hr  r.  . 


(2.2) 


vmere  m is  the  rest  mass  of  the  electron,  M that  of  the  atomic  nu- 

(19) 

C.L.3US.  A usual  approximation'  " to  H^  consists  in  neglecting  the 

off-diagonal  terms  _ II  / V7  .T7  in  Eq.  (2.1).  V(e  have 

2H 
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in  a.u.  (2,3) 


If  the  energy  units  are  now  reduced  "by  a multiplicative  factor  of 

(1_^)  , v/e  obtain  the  modified  atomic  units  (m.a.u.)  for  the  nu- 
2M 

cleus  in  consideration'  . The  new  Hamiltonian  H is  given  hy 


and  we  have 

H'  (a.u.)  = H (m.a.u.) 


(2.4) 


We  shall  use  H throughout  this  calculation. 


2.2.  Basic  Slater  Determinants  and  Spin-orhitals 

(21) 

We  consider  only  normalized  Slater  determinants  D'  ' 

“U_^(_i)?p  -,A=A  = 

III 

made  up  of  orthonormal,  symmetry  adapted  spin-orhitals  ^(^)  in 
v;hich  jjl  = y^(i  ,^,m^  ,m^)  : 


<2-5) 


The  set  [ Y.  j consists  of  normalized  spherical  harmonics  in  the 
Dirac  phase  convention''^^  , and  \^^|are  the  usual  spin  functions 
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cx 


, j3  with  eigenvalues  . The  set  of  functions  is  expanded 

n terns  of  norm 

j = 1,2,. ...Jd) 


( 22 ) 

in  terns  of  nornalized  Slater  type  orhitals  (STO's)  ' . with 


J(^) 

“ie  “ 23 


(2.6) 


(2.7) 


(2Z  . )(2^"^2nj  +3)  /(2E+2n.  +2) ! 
J£  J 


and  v/e  have  further  the  normalization  condition 


/R.,(r)  R.,|j(r)  r2  dr  = 


The  n.*s  are  integers  and  n .'^0  ; the  values  n.  = 0,  and  n.  = 1 are  alv/ays 

J J J J 

included.  Higher  terms  are  introduced  v/hen  the  number  of  orbital  ex- 
ponents for  each  value  of  t is  less  than  three.  There  is  no  perma- 
nent reason  for  this:  it  arises  from  previous  experience  of  other  au- 
thors^^^^  as  well  as  our  own.  Anyiriray,  there  is  not  much  point  in 
simplifying  further  the  expansion  (2.6)  at  this  stage;  v;e  shall  come 
back  to  this  point  later,  in  Chapter  III.  The  powers  of  r in  Eq. 

(2.7)  are  » the  2 being  present  in  order  to  make  "be  in 

the  domain  of  the  Hamiltonian:  specifically,  so  that  some  integrals 
involving  the  operator  l/r^g  '^P*  starting  orbital 

exponents  Z^^  are  selected  by  trial  and  error.  In  general,  v/hen 
optimization  (w'ith  respect  to  total  energy)  is  desired,  a brute 
force  method  is  employed.  The  previous  experience  of  other  workers 
indicates  that  a careful  optimization  does  not  pay^^^^  , and  that 
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once  some  reasonalDle  Z are  obtained,  it  is  much  better  to  intro- 
duce  additional  configurations.  However,  this  is  true  up  to  a certain 
stage,  after  v/hich  the  Cl  expansion  appears  to  converge  very  slowly. 

Then,  it  is  desirable  to  think  about  an  orbital  exponent  optimization, 
which  does  not  need  to  be  done  by  brute  force  alone,  since  by  now, 
much  more  is  known  about  the  regions  where  the  electronic  density  prob- 
ability is  localized.  Further  comments  shall  be  made  in  Chapter  IV. 

In  order  to  get  started  on  our  calculation,  and  because  an  initial 
optimization  of  orbital  exponents  was  impractical  with  an  IBM  709,  we 
v;ere  compelled  to  assume  that  orbital  exponents  which  give  a good 
HF  ground  state  could  also  be  used  to  form  an  STO  basis  which  is 
essentially  correct.  For  orbital  exponents  corresponding  to  orbitals 
with  an  2.  value  not  represented  in  the  ground  configuration,  an 
"excited  configuration"  was  considered,  in  which  some  ground  state 
orbitals  have  been  replaced  by  those  in  question,  and  the  energy  vas 
optimized  while  keeping  the  remaining  ground  state  orbitals  fixed: 
this  vas  done  by  V/atson  in  his  treatment  of  Summarizing,  two 

main  elements  are  employed  in  the  preliminary  search  for  the  one  electron 
functions:  the  HF  orbitals,  and  an  associated  STO  basis  defined  in 
terms  of  the  orbital  exponents  and  of  the  n . power  exponents.  The 
problem  consists  then,  in  obtaining  an  orthogonal  basis  of  a dimension 

smaller  than  that  of  the  STO  basis  and  which  is  convenient  for  a 

(23) 

Cl  treatment.  Watson  used  a method  v/hich  employed  succesive 
Schmidt  orthogonalizations  of  the  STO's,  combined  with  second  order 
rrurbation  theory  estimates  of  the  Cl  effects  introduced  by  some 
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"selected  configurations."  For  more  details,  the  reader  is  referred 
to  Watson’s  paper.  It  is  interesting  to  note  that  additional  func- 
tions of  a given  symmetry  type  were  introduced  in  his  calculation 
until  further  add^vions  contributed  to  the  energy  with  less  than 
0.00005  first  sight,  the  method  may  seem  quite  laborious; 

however,  the  v/hole  process  can  be  carried  out  automatically  on  an 
electronic  computer,  in  a very  simple  way. 

An  apparent  shortcoming  arises  in  the  choice  of  the  "select- 
ed configurations"  to  be  used  in  each  case.  It  presupposes  a know- 
ledge of  which  configurations  are  dominant,  which  may  or  may  not  be 
correct.  As  a result,  the  "selected  configurations"  do  become  artifi- 
cially dominant.  By  including  a reasonable  variety  of  "selected  con- 
figurations" this  difficulty  is  overcome.  Second  order  energy  esti- 
mates are  highly  desirably,  as  they  are  a cheap  alternative  to  a rig- 
orous total  energy  criterion,  which,  after  all,  does  not  necessarily 
lead  to  a better  choice  of  one  electron  functions. 

izh) 

However,  from  the  results  of  V/eiss’  subsequent  work  on  Be 

(20 

and  from  the  recent  investigation  of  Killer  and  Ruedenberg  , v/e 
developed  the  conviction  that  V/atson's  one  electron  functions  could 
be  significantly  improved  simply  by  modifying  the  expansion  coeffi- 
cients while  leaving  the  STO  basis  intact.  It  seemed  possible 

that  some  important  configurations  be  missing  in  Watson's  wave  func- 
tion. Obviously,  v;e  needed  to  find  a better  method  for  the  search  of 
spin-orbitals. 

Our  conviction  had  less  to  do  with  a scientific  analysis 
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than  with  a genuine  interest  in  the  general  applicability  of  the  C.I. 

method.  Weiss'  basis  was  not  orthogonal  and  he  reached  a limit  where 

the  mathematical  stability  of  his  secular  equation  made  it  impossible 

to  add  generously  some  much  needed  configurations.  The  transformation 

to  an  orthogonal  spin-orbital  basis  would  have  increased  the  number  of 

(Z6) 

configurations  in  a manner  probably  unnecessary  and  impractical'  . 
Thus,  as  Weiss  stated  in  his  paper,  his  results  could  not  be  improved 
significantly  along  the  same  line  of  work, and  he  drew  a pessimistic 
conclusion  with  which  we  could  not  agree.  On  the  other  hand, Miller 
and  Ruedenberg's  approach  seems  to  be  general  in  theory  as  well  as 
in  practice,  and  many  stimulating  results  have  been  obtained  and 
are  to  be  expected  in  the  future.  Their  wave  function  is  just  an 
intelligently  constrained  C.I.  wave  function.  Since  there  is  always 
time  to  decide  at  which  stage  and  how  to  impose  constraints  on  a 
C.I.  wave  function,  we  choose  to  follow  the  lines  of  the  conventional 
C.I.  method. 

2.3.  Basic  functions 

A pure  basis  is  one  whose  elements  are  eigenfunctions  to 
some  set  of  quantum  mechanical  operators,  ^flien  these  operators  com- 
mute with  some  Hamiltonian,  the  introduction  of  symmetric  projection 

operators  brings  about,  not  only  a method  for  the  obtaining  of  the 
(27) 

pure  basis  ' , but  also  some  important  simplifications  in  the  secu- 
lar equation  and  in  the  matrix  elements  corresponding  to  that  Hamil- 
tonian^^^^;  this  is  the  reason  for  using  a pure  basis  (the  pure  basis 
alone  simplifies  the  secular  equation;  when  it  is  obtained  by  means 


I 
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of  symmetric  projection  operators,  the  matrix  elements  are  also  sim- 
plified). 

In  the  atomic  case,  the  electrostatic  Hamiltonian  H commutes 

with  L^,  L , S^,  S (squares  and  z-components  of  the  total  orbital  and 
z z 

spin  angular  momentum)  and  v/ith  P , the  parity  operator.  Thus,  we 
use  a basis  of  eigenfunctions  to  (l^,  L , S^,  S , r ] . When  looking 
for  the  variational  solutions  to  H corresponding  to  a given  set  of 
eigenvalues  to  [l^,  L , S^,  S , 'P  ] , we  need  to  concentrate  only 
on  the  subsecular  equation  (of  infinite  order,  hov;ever)  associated 
with  the  above  mentioned  eigenvalues. 

The  projection  operators  0 considered  here,  were  introduced 
by  Lowdin^^^^;  they  are  idempotent,  hermitian  and  symmetric  (so  that 
they  commute  with  H): 

0^  =0  = 0^  ; (3(x^,X2,...,x^j)  = P(9(x^,X2,...,Xj^) 

and  may  be  expressed  as 

0^  =Tf  (2.8) 

where  A fulfills 

A 0j  = 0^ 


and 


(2.9) 


that  is,  iD^  selects  the  component  of  0 which  is  an  eigenfunction  to 
A with  eigenvalue  We  will  be  concerned  only  with  A=  and 
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A=  S2  , since  v;hen  symmetry  adapted  spin-orbitals  are  used,  the 

determinantal  fiinctions  are  already  eigenfunctions  to  L , S and 

z z 

parity.  The  projection  operator  employed  will  he  the  product  of 
(0'-^  times 0 2 


l2 


<3l.s  ' Os 

The  functions  ^ upon  which  the  projection  operator  acts,  are 
a set  of  orthonormal  Slater  determinants  D ; K denotes  the  particu- 
lar  configuration  and  o(  classifies  the  possible  determinants  in  con- 
figuration K (a  configuration  is  defined  by  the  ordered  set  of  all 


(i,E;  numbers  associated  with  the  spin-orbitals  which  form 

a given  determinantal  function). 


). 


Some  configurations  allow  for  only  one  linearly  independent 

projected  function  corresponding  to  a given  set  of  L,L  ,S,S  values 

z z 

f \ 

(each  configuration  has  already  a given  parity  equal  to  ^ ' 

Such  pure  functions  are  called  non-degenerate,  '.'/hen  more  than  one 
linearly  independent  function  can  be  obtained,  these  functions  are 
called  degenerate. 

The  orthonormal  set  of  pure  functions  can  be  expressed 


as 


S^  r T) 

§(p)  = 0 

^L,S 


Y\U  Ko< 


3 c^P^ 


(2.10) 


In  addition,  the  following  relations  are  fulfilled 

2 


O 


L,S 


= 0 


L^,S^ 

L,S 
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The  superscripts  p and  q denote  the  p^^  and  q^^  elements  of  the  ortho- 


normal pure  functions  associated  with  the  configuration  K and  L and 


the  above  equations  is  to  obtain  neat  expressions  for  the  matrix 
elements  of  H.  We  have 


in  that  the  summation  starts  with  ^=1  instead  of  p=q  , the  definitions 
of  the  b’s  and  c’s  being  somewhat  different.  Our  equation  contem- 
pltes  a simpler  computer  program.  When  working  by  hand,  or  in  the 
case  of  large  degeneracies,  the  matrix  elements  are  more  readily  com- 
puted by  means  of  the  Lowdin’s  formula. 

In  addition,  when  gathering  data  to  be  fed  into  the  computer, 
it  is  worthwhile  to  present  them  as  simply  as  possible,  to  avoid 
mistakes  which  otherwise  may  remain  undetected.  Thus,  our  starting 
set  of  pure  functions  shall  be,  instead: 


(2.11) 


where  we  have  made  use  of  the  relations 


(0^  =0=0^  and  (9h=H(^. 

Equation  (2.11)  differs  from  the  one  given  by  Lowdin^^"^^, 


(2.12) 
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where  is  not  necessarily  a normalization  or  a phase  factor.  In 
this  way,  we  leave  the  questions  of  sign  and  normalization  to  the 
computer,  v/hich  transforms  the  into  and 

as  follows;  the  set  is  obtained  from  ^>7  a Schmidt 

orthogonalization 


$(1)  = S<1> 

$(")  = 4"f.g  ($<P)| 

P 


(2.13m) 


The  Eq.  (2,13m)  may  now  he  written  more  explicitly  as 


^(m)  ^ ^(m)f  ^ ^(p)ij^  d^"'H)(f^D  c^P^) 

^K«c°Kcc  1^  K(i%  ^ 


K 


‘f:  4^’4"’>4v’  *4"’i 

K(X  Kot  Kj  J 


“(ra) 


The  unnormalized  coefficients  c^^  are  given  by 


""Kj}  °K«  V ^ 


,(m) 


and  the  normalized  coefficients  c4“‘'  are  then 

H 


(2.1^) 


4:’ = 4:’ 


(2.15) 
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To  obtain  the  h 


(m) 

Kif 


coefficients,  equation  (2.13m) 


is  rewritten  as 


Km 


K Km 


(2.16) 


Thus,  we  have  the  unnormalized  coefficients  h 


>;(ra) 


K> 


for  K<ni 


lim)  ^ j^(m) 
Km  K 


and.  the  normalized  coefficients  are  then 


= t«(g{4“>}b-4 


The 


constant  is  easily  obtained  from  Eq.  (2.12)  as 


(2.17) 


Ml 


(m) 


'n. 


=Ei4:'i^/4:’ 


p(=i 


(2.18) 


2.4.  Basic  Confi^irations 

In  the  previous  section,  the  basic  many  electron  functions 
were  introduced  in  a most  general  way.  It  is  our  purpose  now  to 
deal  with  the  question  of  degeneracies  and  of  how  to  obtain,  in 

]• 


practice,  the  coefficients  ^d^ 


(P) 


Koc 


Zh 


2 2 

Today,  computer  programs  are  availalDle  to  carry  out  L and  S 

projections  of  Slater  determinants  made  up  of  symmetry  adapted  spin- 

orditals^^^^  (the  orbital  quantum  number  I must  not  exceed  the  value 

-d  =^,  in  case  a machine  of  reduced  core  storage,  like  the  IBM  709 » 

is  being  used;  some  special  cases  involving  higher  t values  are  also 

feasable  in  that  machine).  Furthermore,  the  degeneracies  of  the  terms 

that  appear  in  each  case  have  been  well  studied^ and  they  can  be 

evaluated  by  means  of  existing  computer  programs^^^^.  In  spite  of 

1 

this,  we  thought  that  our  application  to  even  S states  of  Be  was 
so  simple  that  it  would  be  far  more  amusing  to  do  this  by  hand,  in- 
stead of  spending  our  time  trying  to  make  operational,  computer  pro-- 
grams  v/hich  had  been  written  by  others.  It  is  instructive  to  do  this 
once  in  a life  time;  next  time,  we  hope  to  use  the  electronic  com- 
puter. 

The  determination  of  degeneracies  is  very  simple.  First,  all 
equivalent  spin-orbitals  (those  with  the  same  (i,E  ) numbers)  are 
grouped,  and  the  possible  terms  of  each  group  are  found  according  to 
the  Pauli  principle:  no  more  than  one  electron  for  each  spin-orbital. 
Then,  all  possible  terms  arising  from  all  groups  of  equivalent  elec- 
trons are  coupled  together,  and  without  restrictions,  according  to 
the  rules  of  addition  of  angular  momenta.  For  instance,  the  configu- 
ration  (s)  (p)  has  two  groups  of  equivalent  electrons.  The  arrange- 
ment  (s)'^  can  give  only  a S term,  while  (p)  gives  three  terms;  D, 

^ and  ^S.  Of  the  latter,  only  the  can  be  coupled  with  the  other 
1 1 

S to  give  a four  electron  S.  Thus,  the  degeneracy  is  equal  to  one. 
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The  configuration  ss'pp'  , instead,  is  douhly  degenerate:  s and  s' 

may  he  combined  to  give  two  terms:  and  while  from  p and  p' 

one  may  obtain:  ^S,  ^P,  and  h.  The  two  possible  four 

1 1 

electron  terms  arise  from  combining  (ss*)  ^ with  (pp')  ^ and 

3 3 

(ss* ) ^ with  (pp' ) 

As  for  the  angular  momentum  and  spin  projection,  Lowdin 

expands  in  terms  of  L =L  +iL  and  L =L  -iL  (step-up  and  step- 

+ X y - X y 

dovm  operators)  and  similarly  for  S^: 


2 2 
= L _ L 

- + z z 


2 2 
S = S - s 

- + z z 


(2.19) 


By  substituting  from  Eq.  (2.19)  into  Eq.  (2.8)  he  obtains,  in  gene- 

(31) 


ral 


= (2M+1)  (-1)  -i2= ^ ^ 

^'‘‘»^'^z  (M-M  )!  -P  =0  ^>!(2M+9+l) 


(2.20) 


where  M stands  for  either  L or  S and  the  caret  distinguishes  between 
operator  and  quantum  number;  Mraax  depends  on  the  particular  configu- 
ration and  is  obtained  according  to  the  rules  of  addition  of  the 
general  angular  momenta.  Eor  the  particular  case  of  the  terms, 
we  have 


''0  A 0 

M- 


'»>!  (>)+!)! 


(2.21) 


In  Table  1,  we  give  the  non-degenerate  configurations  employed 
in  this  work,  together  with  the  projected  functions  and  the  d coeffi- 
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cients.  Table  2 shows  the  doubly  degenerate  configurations  and 
Table  3,  the  triply  degenerate  ones.  In  the  columns  labeled  d'  ' y 
d(2)^  ^{3)  corresponding  to  different  functions  of  the  degenerate 
set,  we  also  have  indicated  which  Slater  determinant  was  projected 
in  each  case.  ¥e  alv/ays  performed  the  L projection  first,  and  then, 
the  spin  projection.  Those  configurations  having  only  s-type  orbitals 
need  a spin  projection  only,  and  they  were  projected  exactly  accor- 
ding to  Eq.  (2.9),  i.e.  =1  (see  Eq.  (2.12)).  The  latter  results 
were  in  turn  applied  directly  to  the  projections  of  the  more  gene- 
ral configurations.  In  the  normalized  Slater  determinant 
the  letters  i,j,k,l  stand  for  the quantum  number,  the  superscripts 
a,b,c,d,  for  m ; the  subscripts  w,x,y,z  stand  for  the  quantum  num- 
ber  i:  in  these  tables  they  are  included  only  to  distinguish  between 
groups  of  equivalent  electrons  and  they  are  arbitrary  otherwise; 
finally,  the  orbitals  to  the  left  of  the  bar  are  associated  with  o< 
spin,  those  to  the  right,  with  a spin. 
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TABLE  1 


NON  deg©;erate  configurations 


Configuration 


Slater  determinants 


(s  s Is  s ) 

W X ' W X 


(s  )^s  s 
w X y 


( s s I s s ) 
w x'  w y 

( S S I s s ) 
w y'  w X 


(s)^(p)^ 


(s  P°|S  P°) 
(s  p^|s  p^) 
(s  /|S  p^) 


s s (p)‘ 

W X ^ 


(s„P^|SiP'') 

(sy  isy> 


is  ) p„Pjt 


(s  p2is  p2) 


ry  1-  i'x' 

,0  o> 


(s  Pxl=  Pw> 

(S  pjis  p'J) 

(s  p'^is  pi) 
(s  pjis  pj) 


d coefficients 


1 


1 

-1 

-1 


1 

1 

-1 

-1 


-1 


-1 


1 

1 

-1 


-1 


-1 


HlM  HlW 
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TABLE  1 continued 

Configuration 

Slater  determinants  d coefficients 

( s 1 s 

X W \ 

(p)** 

(pVip¥)  1 

(pS'^IP^p®)  1 

(pS'^IpV^) 

(p  )3p 
v;  X 

(plp^lpOp’!)  1 

W w'  W X 

(pop'^lp^p®)  1 

W X w w 

(p^p’^lp^p"^)  -1 

W W W X 

(plp'l  iplp'^)  -1 

VJ  X w w 

(pj,p°lpyp°) 

W W W X 

(p"^p°lp^p°)  -1 

W X w w 

(Pypylp^)  ^ 

(p^p'^lp^p*^)  1 

W X w w 

(pOpl|pV) 

W w'  W X 

(pVipS*^) 

W X w w 

(p^p'^lp'^p^)  1 

W W W X 

(P'^P^IPV)  ^ 

W X w w 

(s)2(d)2 

(s  d^is  d®)  1 

( s d^  1 s d'^ ) -1 

(s  d'^ls  d^)  -1 
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Configuration 


s(p)^d 


TABLE  1 continued ' 

Slater  determinants  d coefficients 


(s  d^is  d"^)  1 

( s d^l s d^)  1 


(sd^isdO) 

(s  d'J^ls  dj^) 
(s  d’^ls  db 
(s  d^is  dj,) 

(sdjiedt) 

(s  d^ls  d,,) 
( s il  IS  dj) 
(s  d^ls  d'^) 


/ 0 OjO\ 

(s  p |p  d ) 

(p  d IS  p ) 

f 1 1/2  X 

(s  p |p  d ) 

t 1/2.  1^ 

(p  d 1 s p ) 

/ -1  -l-^2x 

(s  p IP  d ) 

/*1^2.  >lx 
(p  d 1 s p ) 

(s  p |p  d ; 

f 1/1. 

(p  d 1 s p ) 


4 

4 

m 

m 

ri4 

-fTz 

-(12 
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Configuration 


S 


TABLE  1 continued 


Slater  determinants  d coefficients 


(s 

p°dl) 

s i) 

-/IT 

(s  pi 

p°dl) 

-{iz 

iU 

s p^) 

(s  p^ 

p^d^) 

.{Iz 

(pldl 

s p°) 

-(IT 

(s  pi 

p^dO) 

2 

(pld° 

s pM 

2 

(s  pi 

p'dO) 

2 

(pldO 

s p^ ) 

2 

(s  d°  d°d°) 

W W X 

2 

(d°d°  s d°) 
W X w 

2 

(s  d^  d^d°) 

W W X 

-1 

( d^  d*^  s d^  ) 
W X w 

-1 

(s  d^  d^d°) 

W W X 

-1 

(d^d°  s d^  ) 
W X w 

-1 

( s d°  d^  ?■ ) 

W W X 

-1 

(d^d^  s d°  ) 
W X w 

-1 

( s d?'  ^ d^  ) 

W W X 

-1 

d^  s d^  ) 
W X w 

-1 

(s  ^ (f  d^  ) 

--  -1 

w w X 
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TABLE  1 continued 


Slater  determinants 


d coefficients 


(dVis  d^) 

W X w 

-1 

(s  d°|d‘%^) 

W W X 

-1 

(dV|s  d°) 
W X w 

-1 

(s  d^ld^d^) 

W W X 

rr 

(d^d^is  (f^) 
W X w 

c 

(s  d^id^d^) 

W W X 

(d^d^fs  d^) 
w x'  w 

(s  d^ld^d’^) 

W W X 

(d^d‘^|s  d'^) 
w x'  w 

nIT 

(s  d^ld^d'^) 
vr  V/  X 

(T 

(d^d'^ls  d^) 
w x'  w 

r 

(s  d^id^d*^) 

V/*  W X 

(d^d'^is  d^) 
W x'  V/ 

(s  d'^ld'V) 

W W X 

(d'Vis  d'^) 
w x‘  w 

-2 

(d^a°|S  d®) 

-2 

(s  d^ld'^d^) 

W W X 

-2 

(d'^d°|s  d^) 
w x'  w 

-2 

(s  d°|d^d^) 

W W X 

-2 

(d^d^ 
w X 


s d®) 
w 


-2 
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Configuration 


TABLE  1 continued 


Slater  determinants  d coefficients 


(s  d2|d°d"2) 

W W X 

-2 

(d°d"^|s  d^) 
w x'  w 

-2 

(s  d'^id^d^) 

w'  W X 

-2 

(dVis  d’2) 

W X w 

-2 

(s  d°|d^d^) 

W W X 

-2 

(d^^d^ls  d°) 
w x'  w 

-2 
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TABLE  2 


DOUBLY  DEGEIIERATE  CONFIGURATIONS 


Configuration 

Slater  determinants 

(Sy^SylS^S^) 

(s^SxISyS^) 

^VyJ^x^z) 

1/3 

1/6 

(s^SxISySj 

1/6 

1/3 

(s^szis^sy) 

1/3 

1/6 

(SyS.lS^Sx) 

1/6 

1/3 

(S^SglSjjSy.) 

1/6 

-1/6 

1/6 

-1/6 

(swP?'SxP°) 

(SwP?|SxP?^ 

s s p p 
w X y z 

(s^I^ls^pO) 

1 

1/2 

V 

(SwP^SxpJ) 

1/2 

1 

(s^pjis^pl) 

-1 

-1/2 

(s^pjis^l^) 

-1 

-1/2 

-1 

-1/2 

(SxP^IS^^) 

-1 

-1/2 

(Sw^xipyi?) 

1/2 

-1/2 

(lyPz  1^?^) 

1/2 

-1/2 

1%^  ^ 

-1/2 

1/2 

^%Pz  ^ 

-1/2 

1/2 

^ l^Pz  ^ 

-1/2 

1/2 

^^yPz  ^ 

-1/2 

1/2 
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TABLE  2 continued 

Slater  determinants  d^^^  d^^^ 


/ 0. 

(s  p IS  p ) 
X z w y 

1 

1/2 

(s  p IS  p ) 
X y w z 

1/2 

1 

(s  p is  p ; 
w z x^y 

-1/2 

-1 

/ -1  Iv 

(s  p IS  p ) 
x^y*  wz 

-1/2 

-1 

/ 4 In 

(s  p |s  p ) 
w z ' X y 

-1/2 

-1 

/ 1 4 \ 

(s  p )S  p ) 
x^^y'  w z 

-1/2 

-1 
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TABLE  3 


TRIPLY  DEGENERATE  CONFIGURATIONS 

d(l) 

Configuration  Slater  determinants  ^ ^ , 

'❖»y>  <P»^y> 


‘P»P:<IP»Py’ 

7 

-3 

1 

(pV|p¥) 

wx'wy 

-3 

12 

0 

(p^p°ip°p^) 

w w-  X y 

5 

0 

2 

J 

(P*^P^|P^P°) 
w y'  w X 

7 

-3 

1 

(P°P^|P*^P°) 

wx'wy 

7 

-3 

1 

(p^p^ip*^p^ ) 

^w  yi^w^x 

7 

-3 

1 

(p*^P^  1P^P°) 
w X ' w y 

7 

-3 

1 

(p^p*^  IP^P^ ) 
^w^yl^w^x 

7 

-3 

1 

/-l^O.  0 lx 
^rf^x  ^w^y 

7 

-3 

1 

(p*^p^  iP^P*^ ) 
^w^y  l^w^x 

7 

-3 

1 

,0-1  1 On 

(PxPylPwPw) 

5 

0 

2 

,0-1  1 0 V 

^PwPw'PxPy^ 

5 

0 

2 

,10  0-lv 

^PxPy'PwPw^ 

5 

0 

2 

,10  -1  On 

(PwPw'Px^^^^ 

-5 

0 

-2 

(p^P^  IP^P^  ) 

-5 

0 

-2 

(p^p^  iP^p^ ) 

^w^w ' x^y 

-5 

0 

-2 

(P°P^  IP^P^ ) 

-5 

0 

.2 

(p^P^  |P°P^ ) 
w y ' w X 

2 

-3 

-1 

(p^p^  IP^P^ ) 
'^w^y 

2 

-3 

-1 
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Configuration 


TABLE  3 continued 


Slater  determinants 


2 

-3 

-1 

2 

-3 

.1 

2 

-3 

-1 

2 

-3 

-1 

2 

-3 

-1 

2 

-3 

-1 

(p2;>2pi> 

-3 

12 

0 

/_0  "1  0 1 X 

<VxlVy> 

-3 

12 

0 

,01  0-lx 
(P„Py|P„Pj) 

-3 

12 

0 

,10  -10, 
<VxlR,Py) 

-3 

12 

0 

,-10  10, 
'VylPwPx> 

-3 

12 

' 0 

,-l  0 10, 

-3 

12 

0 

(Pip“,iip“) 

-3 

12 

0 

<pipxipi^y) 

-4 

-9 

-1 

(piPylPjp^) 

-4 

-9 

-1 

(P^P^IP^P^) 

-4 

-9 

-1 

-4 

-9 

-1 

<pipilPx^^ 

-5 

0 

-2 

(ppjipi^i) 

-5 

0 

-2 

(ivpii^pj’ 

5 

0 

2 

(p^p^  IP^P^ ) 

^x^y'^w^w 

5 

0 

2 

(pipJ4j^> 

1 

_9  -- 

1 
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Configuration 


TABLE  3 continued 


Slater  determinants 


,(1) 


,(2) 


,(3) 


( 1 ) 
w X w y 

1 

-9 

1 

(p^P^  IP^P^ ) 
w y ' w X 

1 

-9 

1 

(P^P^  IP^P^ ) 
w X ' w y 

1 

-9 

1 

(P^P^  IP^P^ ) 
w y ' w X 

-9 

6 

0 

-9 

6 

0 

<5i4litPy> 

-9 

6 

0 

<P>y'P>x> 

-9 

6 

0 

/ 0 0 0 Ox 

(p„Pxip„iV) 

-6 

-6 

0 

/ 0 0 0 Ov 
‘Vy'Vx> 

-6 

-6 

0 

(pSp>S4) 

IP^P^ ) 
^w^x  ^vrx 

7 

-3 

1 

(p^p^  iu*^-n^ ) 

^w^x'*w^x 

12 

0 

/ 1„0  0 "1  \ 

5 

0 

2 

/_0“1  _1_0\ 
(p  p IP  p ; 

7 

-3 

1 

/_0  1 “1_0  X 

<Vx'Vx> 

7 

-3 

1 

(p^p^i p^p^ ) 
^w^x'  ^w*^x 

7 

-3 

1 

(p'^p^l  p^p*^) 
^x^x' 

5 

0 

2 

/^O-l.U^Ox 

5 

0 

2 

,10  0-lx 

^VxIVw^ 

5 . 

0 

2 

,1-1  0 0 X 

^Vx'Vx^ 

2 

-3^_ 

-1 
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TABLE  3 continued 


Slater  determinants 

d'l> 

d^2^ 

d<3) 

(pViP^p'^) 

W X W X 

2 

-3 

-1 

(pVip'S^) 

W X W X 

2 

-3 

-1 

(pS^pV) 

W x'  W X 

2 

-3 

-1 

(p'^p^lp^p^) 

-3 

12 

0 

-3 

12 

0 

-3 

12 

0 

<p»yx> 

_4 

-9 

-1 

'pyx'p>x> 

-4 

-9 

-1 

(p^p'^ipV^) 

^w^w' ^x^x 

-5 

0 

-2 

,1-1  1 -1, 
'I’xI’xIVw' 

-5 

0 

-2 

-9 

^ -1  1 *1  1 ^ o L 

(p  p It)  p ) -9  o 

(pSp>S^$) 
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CHAPTER  III 


CALCULATION 
3.1.  The  STO  Basis 

We  employed  the  STO  basis  previously  used  by  Watson^^^^,  to 
expand  s-,  p-  and  d-type  orbitals.  In  his  calculation,  V/atson  consid- 
ered f-  and  g-type  orbitals  also;  these  were  not  considered  here. 

V/atson's  choice  of  orbital  exponents  was  based  on  previous 
experience  from  his  own  calculations  on  the  He  ground  state^^^.  We 
checked  his  HE  orbitals  by  running  Roothaan*s  atomic  SCF  program^ 
on  the  IBM  709  computer,  v;hile  keeping  V/atson' s ST0*s  fixed;  we  found 
them  correct  to  all  ciphers  in  the  energy  and  at  least  to  five  fig- 
ures in  the  orbitals.  The  STO  basis  is  given  in  Table  A,  for  com- 
pleteness. As  a result  of  limiting  ourselves  to  only  two  orbital  expo- 
nents for  each  type  of  orbital  symmetry,  it  was  possible  to  include 
several  powers  of  r without  much  programming  effort.  Watson's  com- 
promise of  two  orbital  exponents  and  n.  values  up  to  n . =k  for  s-type 

J J 

orbitals,  is  very  similar  to  the  option  chosen  by  Weiss,  who  combined 

only  one  n.  value,  either  n =0  or  n =1,  with  each  of  four  orbital 
J J J 

exponents. 
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ho 


TABLE  h 


THE  STO  BASIS 


J=1 

2 

3 

h 

5 

6 

7 

8 

9 

n .=0 
J 

0 

1 

1 

2 

2 

3 

3 

4 

jl 

6.0 

1.0 

6.0 

1.0 

6.0 

1.0 

6.0 

1.0 

6.0 

J2 

9.0 

1.5 

9.0 

1.5 

9.0 

1.5 

9.0 

J3 

12.0 

2.0 

12.0 

2.0 

12.0 

3.2.  The  Search  for  s and  p Orbitals 


According  to  the  scheme  outlined  in  Section  1.3>  we  performed 
Cl  calculations  based  on  Eqs.  (1.9)  and  (1.10)  and  with  s-  and  p-type 
orbitals  only.  Two  different  sets  of  spin-orbitals  were  employed 
as  a basis:  set  A was  made  up  of  the  NSO’s  of  Watson's  wave  function 
^^3)  as  calculated  by  Barnett^^^^;  set  B consisted  of  a trial  basis: 
the  orbitals  s^,  s^  and  p^  were  taken  to  be,  in  order,  the  Is,  2s 
and  p^  of  Watson.  Additional  orbitals,  up  to  s and  p^,  are  cons- 
tructed  by  successive  Schmidt  orthogonalization  of  the  STO  basis 
according  to  the  order  j=2,3,5>6,7  for  s-type  orbitals  and  j=l,3»^»5»6 
for  p-type  orbitals. 

The  results  are  displayed  in  Table  5*  The  energy  differences 
between  Cl  results  obtained  with  the  two  basic  sets  are  not  considered 
to  be  significant.  In  Table  12,  the  occupation  number  spectra  of 
some  of  these  wave  functions  are  shown.  The  resemblances  between 
the  occupation  numbers  of  the  K and  L shell  Cl  NSO's  of  set  B with 
some  of  those  from  the  Weiss  wave  function,  are  striking.  The  K and 
L shell  NSO's  are  very  complementary,  in  the  sense  that,  when  one  of 
them,  say,  from  the  K shell  wave  function,  has  a relatively  high 
occupation  number,  the  corresponding  NSO  for  the  L shell  has  a very 
small  one,  and  vice' ver;sa(  the  NSO's  from  both  K and  L shell  Cl  were 
classified  by  comparing  their  plottings  with  those  from  the  KL  shell 
Cl;  the  occupation  numbers  served  as  a useful  guide). 

Since  the  NSO's  from  the  KL  shell  Cl  are  orthogonal  and  since 
they  match  quite  well  their  K or  L counterparts,  the  new  spin-orbital 
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TABLE  5 


ElIERGY  VALUES  AT  DIFFERENT  STAGES  OF  THE  SEARCH  FOR  p ORBITALS 


Wave  function 

Comments 

Energy 

Watson's  Hartree-Fock 

-14.572986 

K shell  C.I.  with 
Watson's  NSO  "basis 

-14.610596 

Idem  with  trial  basis 
(see  text) 

-14.610624 

*<) 

L shell  C.I.  with 
Watson's  NSO  basis 

-14.617783 

Idem  with  trial  basis 

-14.617391 

(Si)^(S2)^+(S2>^2I(‘l*j)'^ 

KL  shell  C.I.  with 
Watson's  NSO  basis 

Idem  with  trial  basis 

-I4.65I856 

-14.651691 

Full  C.I.  with 

125  C.I.  with 
NSO's  of  KL  shell  C.I. 

-14.657015 

» Sg*  s^, >P2»P2 
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"basis  v;as  constructed  from  them. 

¥e  had  to  decide  which  set  of  KL  shell  NSO's  to  investi^te 
first,  and  we  chose  the  HSO’s  from  set  A "because  it  seemed  more  conven- 
tional to  do  so,  and  there  were  no  other  reasons  to  prefer  either 
choice.  In  a way,  set  B was  more  interesting,  regarding  the  proce- 
dure "by  means  of  which  it  had  "been  o"btained,  "but  it  seemed  somewhat 
more  risky. 

With  this,  we  concluded  our  search  for  s and  p orbitals;  the 

s and  s„  orbitals  were  to  remain  the  original  HP  ones  (thanks  to 
12 

this,  a cluster  expansion  analysis  of  our  final  wave  function  can  be 
compared  with  the  analysis  of  V/at son's  wave  function  made  by  Sinanoglu 
and  the  others  were  chosen  successively  from  the  highest  occu- 
pation number,  down.  One  can  take  advantage  of  this  situation  and 
define,  for  a given  quantum  number  , the  orbital  number  i according 
to  the  position  of  the  orbital  in  Table  12.  When  the  wave  function 
is  good,  the  occupation  number  spectrum  varies  little  with  further 
improvements;  also,  "crossings"  are  very  rare  and,  for  a given  symme- 
try, the  spacing  is  wide  enough  so  that  a unique  mapping  is  possible 
between  the  NSO's  of  some  good  approximation  to  and  the  NSO's  of 
a still  better  approximation  to  ^ . It  was  with  this  expectation  in 
mind,  that  we  partitioned  a set  of  Slater  determinants  into  two  sets 
L and  L’  (see  Eq.  (1.6)). 

3.3.  The  Search  for  d Orbitals 

A full  Cl  treatment  is  one  in  which  all  possible  configurations 
arising  from  the  spin-orbital  basis  considered,  are  taken  into  account 
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(see  Tables  1,  2 and  3)» 

We  selected  four  s-type  orbitals  out  of  the  six  possible  ones 
and  three  p-type  orbitals  out  of  the  corresponding  five.  The  energy 
of  the  full  Cl  expansion  in  a 4s,  3p  basis  (94  configurations;  125 
terms),  was  calculated  (see  Table  5)»  and  the  result  considered  very 
encouraging  when  compared  with  the  sp-limit  attained  by  Weiss  (see 
Chapter  IV),  The  truncation  of  the  125  Cl  expansion  was  tentatively 
made  by  using  a rough  eigenvector  component  criterion  (see  Table  6) 
which  worked  very  well  in  this  case. 

We  considered  the  four  d orbitals  of  Watson,  and  we  obtained  a 
fifth  one  by  a Schmidt  orthogonalization  of  the  STO’s  (see  Eq,(2,7) 
and  Table  4)  with  respect  to  the  starting  basis.  Next,  some  configu- 
rations including  d orbitals  were  considered,  in  which  only  s^,  s^, 

P^»  P2  p^  were  included.  All  the  possible  combinations  of 
the  form  (s)^(d)^,  (s)^d^d^  and  s(p)^d  were  taken  into  account;  only 

three  configurations  of  type  s(d  )^d  were  included  (see  Table  1); 

w X 

in  this  instance,  as  well  as  in  later  experiments,  their  eigenvector 
components  were  smaller  than  10“'5  and  their  energy  contributions  were 
negligible. 

The  resulting  expansion  contained  143  terms:  the  corresponding 
result  is  shown  in  Table  6,  The  natural  orbitals  of  this  wave  function 
were  calculated  (see  Table  12),  and  we  selected  the  three  d-type 
orbitals  with  highest  occupation  number.  Incidentally,  the  s and  p 
NSO's  were  seen  to  be  identical  with  the  original  ones;  in  the  latter 
ones,  the  former  ones  showed  up  with  expansion  coefficients  equal 
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TABLE  6 

ENERGY  VALUES  AT  DIFFERENT  STAGES  OF  THE  SEARCH  FOR  d-ORBITALS 


Truncated  full  (4s,  3p)  125  C. I. 


Wave  function 

Comments 

Energy 

Only  terms  with  eigenvector 
component  larger  than  5 .10"^ 

85  C.I. 

-14.657015 

Only  terms  with  eigenvector 
component  larger  than  10“^ 

72  C.I. 

-14.657012 

Only  terms  with  eigenvector 
component  larger  than  4 .10"^ 

51  C.I. 

-14.656983 

Inclusion  of  d configurations 

Wave  function 

Comments 

Energy 

(See  text) 


143  C.I. 


-14.660139 


to  0.9992  or  larger;  among  the  p-type  NSO’s,  the  corresponding 
coefficients  were  larger  than  0.999963. 

3.^.  Final  Selection  of  Configurations 

In  order  to  speed  up  the  procedures,  we  left  the  s^,  s^  and 
p^  (NSO's  of  the  KL  shell  Cl  wave  function)  out  of  the  picture.  They 
are  included  now.  The  successive  steps  are  best  described  in  Table  8. 
All  orbital  expansion  coefficients  in  terms  of  the  STO  basis  are  given 
in  Table  7.  An  analysis  of  the  final  results  is  given  in  Chapter  4. 

3.5.  Computing  Techniques 

An  IBM  709  computer  was  used  throughout  this  work.  All  pro- 
grams were  written  in  Fortran  II,  using  double  precision  arithmetic. 

The  basic  two  electron  integrals  were  calculated  by  means  of  the 
recurrence  relations  given  by  Roothaan  et.  al.^35)  and  stored  on 
magnetic  tape  for  use  in  future  calculations.  The  generation  of  the 
155  X 155  Hamiltonian  matrix  of  our  final  wave  function  took  twenty 
minutes.  The  eigenvalue  and  eigenvector  were  found  by  using  a sub- 
routine program  written  by  F.  Prosser^^^  which  employs  the  Givens' 
method;  it  is  written  in  single  precision  arithmetic,  but  includes 
selective  rounding  off  at  all  the  critical  stages. 

A computer  program  was  written  which,  under  certain  conditions, 
improves  an  approximate  eigenvalue  and  eigenvector  to  any  desired 
accuracy.  It  consists  of  an  application  of  the  Gauss-Seidel-Lanczos^^^^ 
iteration  method  to  the  partitioning  technique  equations.  The  nodes 
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of  the  hraclceting  function  are  found  hy  means  of  a direct  interpo- 
lation he  tween  the  value  of  the  bracketing  function  for  an  initial 
energy  and  its  value  after  a 45°  construction I The  bracketing 
function  calculated  at  this  interpolation  point  is  usually  smaller 
than  10"9,  but  since  this  value  may  be  in  error,  we  calculate  the 
final  bounds  for  the  lowest  eigenvalue  of  the  finite  Hamiltonian 
na-trix  in  the  following  way:  let  E and  c be  the  approximate  eigen- 
value and  eigenvector  respectively;  the  residual  vector  is 

TQ  = (H  - E)c 

and  if  c is  normalized,  it  may  be  shown  that^^^^ 

E = E + H 11 

The  Lanczos  method  consists  in  minimizing  the  residual  vector  after 
several  applications  of  the  Gauss-Seidel  iterations. 

The  eigenvalues  and  eigenvectors  of  several  Hamiltonian  ma- 
trices calculated  with  Prosser's  subroutine,  were  checked  against  the 
ones  obtained  with  our  own  program;  the  eigenvalues  were  usually 
accurate  to  8 figures,  while  the  eigenvectors  were  correct  in  3 to  7 
figures. 
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CHAPTER  IV 


ANALYSIS  OF  THE  RESULTS 

In  this  chapter,  we  shall  mention  first  the  salient  features 
of  this  calculation.  Then,  we  shall  discuss  how  to  improve  it.  Noth- 
ing will  he  said  about  the  possibilities  of  obtaining  comparable 
accuracies  for  larger  systems,  since  this  would  be  pure  speculation. 

4.1.  Salient  Features 

The  total  energy  obtained  for  Be  ground  state  is  the  lowest 
one  found  so  far  by  means  of  a rigorous  calculation  (see  Table  9; 
the  estimated  non-relativistic  energy  is  -14.66745  m.a.u.^^^^).  The 
improvement  over  Weiss*  energy  value  is  believed  to  be  due  to  the 
addition  of  more  configurations  containing  p-type  orbitals.  If  only 
s-type  configurations  are  considered,  the  best  energr  for  a wave 
function  with  these  characteristics  is  called  the  s-limit  energy. 
V/eiss*  s-limit  energy  is  -14.59110  m.a.u.  (see  Table  9 of  his  paper); 
ours  is  -14.59073  m.a.u.  The  two  energies  are  so  close  to  each  other 
that  it  is  impossible  to  decide  which  s-type  basis  will  be  more  effec- 
tive in  lowering  the  energy  of  the  final  wave  function.  The  sp-limit 
of  Weiss  is  only  -14.65722  m.a.u.,  which  compares  unfavorably  with 
ours:  -l4. 659132  m.a.u.  ’nHien  adding  configurations  containing  d-type 


50 


51 


TABLE  9 

RIGOROUS  CALCULATION  ON  Be  GROUND  STATE 


Year 

Computer 

Wave  function 

Energy 

i960 

IBM  704 

37  C,I,;  uses,  in  part,  a H,F, 
Basis,  and  up  to  g-type  orhitals 

-14,65740^ 

1961 

UNIVAC 

55  C,I,;  uses  non-orthogonal 
basis  and  up  to  f-type  orbitals. 
Optimizes  non-linear  parameters 

-14,66090^ 

1965 

IBM  7074 

Augmented  separated  pair  expansion 
29  terms,  up  to  d-type  orbitals. 
Optimizes  non-linear  parameters 

-14.66179^^ 

1966 

IBM  709 

80  C,I,;  uses  a mixed  basis  of 
H.E,  orbitals  and  approximate 
NSO's;  up  to  d-type  orbitals 

-14.66203*^ 

155  c,i. 

-14,66213*^ 

^Watson,  R.E.,  Phys.  Rev,  119,  1?0  (i960), 

%eiss,  A,¥,,  Phys,  Rev,  122,  1826  (I96I), 

Q 

Miller,  K.J,,  and  Ruedenherg,  K, , preprint.  Institute  for  Atomic  Research, 
Department  of  Chemistry,  Department  of  Physics,  Iowa  State  University, 
Ames,  Iowa, 

^This  calculation. 
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ortitals,  V/eiss'  improvement  is  of  -O.OO3I7  m.a.u.;  ours  is  only  of 
-0. 00300  m.a.u.  The  improvement  obtained  by  Watson  when  adding  the 
d orbitals  is  just  about  the  same.  Since  we  did  not  include  the 
important  configurations  of  the  type  (p)  (d)  , it  is  probable  that 
our  improvement  will  be  slightly  larger  than  theirs,  though  of  the 
same  order  of  magnitude. 

To  compare  our  wave  function  with  others,  it  would  be  necessary 
to  transform  all  of  them,  including  ours,  into  an  NSO  expansion  form. 
Barnett  et.  al.'  ' did  this  for  V/atson’s  and  V/eiss'  wave  functions. 
Since  our  143  Cl  wave  function  (see  the  end  of  Section  3.3  and  Table  6) 
is  almost  in  natural  form,  it  is  possible  to  make  a comparison  of 
our  expansion  coefficients  with  those  from  the  natural  expansions 
given  by  Barnett  et.  al.  They  are  very  similar,  all  phases  agree, 
and  none  of  Weiss'  configurations  is  missing  from  our  wave  function. 

For  example,  the  first  five  coefficients  of  V^eiss'  expansion  are 
0.95327,  -0. 29626,  -0. 036491,  -0,028171  and  -O.OI7063,  while  our 
coefficients  corresponding  to  the  same  configurations  are  0. 95255, 

-0. 29612,  -0.04o428,  -0. 027849  and  -O.OI6523,  respectively.  The 
small  difference  between  -0. 29626  and  -O.29612,  which  correspond 
to  (s^)  (p^)  , may  be  due  to  the  fact  that  our  s^  and  p^  are  almost 
KSO's;  it  may  also  be  a coincidence. 

In  Table  10,  we  present  some  features  of  our  wave  function 
v.'hich  are  of  wide  interest.  First,  the  truncations  made  in  the 
155  Cl  wave  function,  illustrate  the  goodness  of  the  eigenvector 
component  criterion  in  this  respect.  The  second  series  of  trunca- 
tions is  relevant  both  to  the  cluster  expansion  and  to-  the  problem 


TABLE  10 


ENERGIES  FOR  SOME  TRUNCATED  EXPANSIONS 


Truncations  from  155  O.I. 

Wave  function 

Only  terms  with  eigenvector 
component  larger  than  5 .10“^ 

Only  terms  with  eigenvector  l 
component  larger  than  8 ,10“^ 

Only  terms  with  eigenvector 
component  larger  than  10"^ 

Truncations  from  80  C.I. 

Wave  function 

Triple  excitations  left  out 

Triple  and  single  excitations 
left  out 

Quadruple,  triple  and  single 
excitations  left  out 

Only  (t)  ^(u)  ^ terms 


80  C.I.^ 

Energy 

-14.66202? 

68  C.I.° 

-14.661923 

59  C.I.^^ 

-14.661619 

Energy 

68  C.I. 

-14.661772 

63  C.I. 

-14.661228 

53  C.I. 

-14.658293 

36  C.I. 

-14,652744 

^See  Table  11. 

^In  Table  11, terms  6?  and  69  through  79  are  missing. 

^In  Table  11, terms  58,  59»  ^1  3'inl  ^3  througih  80  are  missing. 
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of  the  convergence  of  the  Cl  method  itself. 

We  shall  define  the  correlation  energy  E as 

corr 

^corr  ” E^j^g^j,^(non-relativistic)  - Ej^j,(non-relativistic)  . 

A rigorous  definition  is  given  hy  Lowdin^^^^, 

In  quantitative  estimations  of  correlation  energies  hy  the  many 
electron  theory^^^^  or  the  infinite  order  perturbation  theory^ 
some  assumptions  are  made  about  the  magnitudes  of  several  terms  appear- 
ing in  the  theory.  To  this  author's  knov/ledge,  a direct  check  on  the 
energy  obtained  has  never  been  done.  Either  the  average  value  of  the 
obtained  wave  function,  with  respect  to  the  Hamiltonian,  can  be  cal- 
culated without  serious  approximations  in  the  integrals,  in  which 
case  the  omission  of  an  energy  check  is  unjustified,  or  some  mathe- 
matically correct  energy  bounds  for  this  average  value  are  not  obtain- 
able, in  which  case  the  method  cannot  be  said  to  be  rigorous.  Tuan 
and  Sinanoglu^^^^  have  reported  a Be  ground  state  calculation  in 
which  99.8?  % of  the  correlation  energy  was  accounted  for.  On  the 
basis  of  an  analysis  of  Watson's  wave  function they  assumed 
that  the  energy  contributions  of  the  single  and  triple  excitations 
v/ere  negligible;  we  see,  from  Table  10,  that  they  contribute  very 
little,  although  significantly. 

Our  quadruply  excited  configurations  differ  sli^tly  from  those 
of  V/atson,  but  their  energy  contributions  are  comparable.  It  is  not 
worthwhile  to  make  a cluster  expansion  analysis  of  our  wave  function 
at  this  moment,  since  v/e  hope  to  improve  it  in  the  next  few  months, 
but  after  that,  v;e  believe  it  will  prove  to  be  useful.  - 
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In  Tatle  11,  the  eigenvector  of  our  80  Cl  wave  function  is 
displayed.  V/e  note  the  relatively  hi^  eigenvector  component  of  some 
single  excitations.  These  are  accurate  to  all  figures  shown  in  the 
table;  they  are  written  here  only  for  the  purpose  of  illustration, 
since,  in  order  to  use  them  in  actual  computations,  more  ciphers  are 
needed. 

In  table  12,  some  fairly  representative  occupation  number  spec- 
tra of  the  wave  functions  mentioned  here  are  shown.  The  rise  of  the 
occupation  number  of  p^  in  the  final  wave  function  can  be  rationalized 
by  noting  that  the  inclusion  of  p^  brings  into  the  Cl  expansion, 
several  configurations  containing  that  orbital  and  with  relatively 
high  eigenvector  component.  Of  course,  the  next  thing  to  be  done 
is  to  employ  our  final  HSO’s  as  a new  basis,  including  s^  and  s^; 
we  shall  consider  this  possibility  in  the  next  section, 

V/e  have  already  pointed  out  in  Chapter  III  some  of  the  uses 
we  have  made  of  Table  12. 

In  Table  13»  we  compare  the  radial  density  distribution  of 
our  final  155  Cl  wave  function  with  those  obtained  by  means  of  a HF 
wave  function.  In  order  to  obtain  the  average  radial  density,  the 
p-type  contributions  must  be  multiplied  by  3,  and  the  d-type  ones, 
by  5.  If  one  compares  the  average  radial  density,  one  finds  that  the 
HF  is  larger  in  the  K shell  region  and  smaller  in  the  L shell  region. 
As  expected,  the  p-type  contributions  are  essentially  localized  in 
the  L shell  region. 


TABLE  11 


AN  80  TERMS  WAVE  FUNCTION  FOR  THE  GROUIID  STATE  OF  THE  Be  ATOM 


Configuration^ 

) 

Eigenvector 

(s^)^S2)^ 

0.952666 

-0.294953 

(s^)^(s3)^ 

-0.039464 

(s^)^s^s^ 

* 

-0.038081 

(S2)^P2)^ 

-0.027936 

(S2)^s^)^ 

-0.017576 

(S2^)^(d^)^ 

-0.016772 

'^2  ^3^4 

0.013383 

2 

(Si)  S2S4 

* 

-O.OO8534 

(Pl)^(P2^^ 

(2) 

-0.008491 

2 

(S2)  P1P2 

-0.008306 

S1S2P1P2 

(1) 

0.008203 

(sJ 

-0.007189 

2 2 
($2)  (^2) 

-0.006559 

Si S2P2P3 

(2) 

0.006490 

S1S2P2P4 

(2) 

0.006082 

-0.005979 

S1S2P1P3 

(1) 

0.005950 

2 

SlS2(Pi) 

0.005668 

(S4)^(pj)^ 

0.005590 

TABLE  11  continued 


Configuration 

Eigenvector 

(33)2(33)2 

-0.005553 

S1S2P1P2 

(2) 

-0.005043 

(s3)2(P3)2 

-0.004801 

SiSgPiPii 

(1) 

0.004455 

S3S^(p^)2 

-0,004131 

(1) 

0.004062 

h^z^3^k 

(2) 

-0.003636 

(2) 

0.003266 

-0.003162 

$132(83) 

-0.003072 

(s^)^(s^)^ 

-0.002836 

Si$2(P2)^ 

-0.002740 

{5^)  8383 

■«■ 

0.002482 

S1S3P1P2 

(1)  *** 

0.002316 

(P^)^(P2;)^ 

(3) 

0.002291 

(1) 

0.002271 

S2S4(Pi)^ 

0.002178 

S182P2P3 

(1) 

0.002131 

(82)^(d3)^ 

-0.002044 

(Pl)\ 

0.001946 

(S3)^pi)^ 

0.001927 

(S;^)^(P4)^ 

-0. 001872 

TABLE  11  continued 


nfiguration 

Eigenvector 

-0.001793 

, ,2 

<Sl)  P3Pi, 

-0.001787 

SjSjCp^) 

O.OOI6A7 

p p 

O.OOI608 

(Pl)  (P3) 

p p 

(2) 

-0.001593 

(S2)^(p^)^ 

-0.001576 

(2) 

0.001544 

=2=3<Pl> 

■«-!«{■ 

-0.001464 

Vkh^z 

(1) 

-0.001355 

2 

(sg)  P2P4 

0.001180 

(1)  •tHHi’ 

-0. 001171 

<=2>  =3=5 
, ,2  .a 

-0.001145 

(s.)  (p„) 

p p 

0.001136 

(S3)  (p^) 

-0.001124 

SjS  (P 

p 

0.001094 

3^33(3  ) 

p ^ p 

-0.001056 

(s^)  (P3) 

-0.001049 

(1) 

0.001037 

/ x2,  ,2 

(sJ  (p^) 
^ 2 

0.001033 

(Sj)  s^s 

0.001020 

2 

Si(pj)  dj 

O.OOO877 

(s  )^s^)" 

p p 

0. 000862 

(S3)  (s^) 

-0.000848 

TABLE  11  continued 
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Configuration 

Eigenvector 

^1^2^3^5 

(1) 

-O.OOO818 

^2  ^ ^3  ^ 

■IHHS- 

-0.000793 

s^s^Cp^)^ 

0.000791 

(S2) 

-0.000762 

^2^3^ 

0.000730 

^1^2^4^5 

(1) 

0.000716 

s^s^p^p^i 

(1) 

0.000675 

S1S5P1P2 

0 

(2) 

0.000609 

c 

(S^)  P;l^P2 

j, 

0.000602 

(Pl) 

0,000599 

(Si)  S2S5 

0.000586 

2 

^1^^^^  ^3 

0.0005^9 

( s ) p p, 

-0. 000541 

^^22^5  2 
(Sj)  {s  ) 

(2) 

0.000528 
-0. 000524 

a 

The  energy  corresponding  to  this  wave  function  is  E=-l4. 662027  a.u. 
h 

The  number  in  parenthesis  indicates  which  term  of  the  degenerate 
configuration  is  used.  The  asterisk  (*)  denotes  a singly  excited 
configuration;  three  asterisks, (***) , a triple  excitation. 


OCCUPATION  NUMBER  SPECTRA  OF  SOME  GROUND  STATE  Be  WAVE  FUNCTIONS 
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TABLE  13 

RADIAL  DENSITY  PROBABILITY  DISTRIBUTION 


Hartree-Eock^  Final  wave  function 


r(a.u.) 

s-type 

p-type 

d-type 

0 

0 

0 

0 

0 

0.01 

0.02053 

0.02055 

0.00000 

0.00000 

0.02 

0.07582 

0.07587 

0.00000 

0.00000 

0.03 

0.1575^ 

0.15756 

0.00000 

0.00000 

0.04 

0.25869 

0.25862 

0.00000 

0.00000 

0.05 

0.373^3 

0.37319 

0.00001 

0.00000 

0,06 

0.49644 

0.00001 

0.00000 

C.07 

0.62437 

0.00002 

0.00000 

0.08 

0.75372 

0.00003 

0.00000 

0.09 

0.88185 

0.00004 

0.00000 

0.10 

1.00818 

1.00667 

0.00006 

0.00000 

0.12 

1.24016 

0.00010 

0.00000 

0.14 

1.44520 

0.00015 

0.00000 

0.16 

1.61726 

0.00022 

0.00000 

0.18 

1.75^86 

0.00030 

0.00000 

0.20 

1.86205 

1.85859 

0.00038 

0.00001 

0.25 

1.98872 

1.98449 

0.00062 

0.00001 

0.30 

1.96401 

1.95911 

0.00085 

0.00002 

0.35 

1.83867 

1.83344 

0.00106 

0.00002 

0.40 

1.65609 

1.65106 

0.00123 

0.00002 

0.45 

1.44901 

1.44468 

0.00138 

0.00002 

0.50 

1.24002 

1.23675 

0.00153 

0.00002 

0.60 

0.86751 

0.86646 

0.00189 

0.00001 

0.70 

0.58938 

0.58945 

0.00244 

0.00001 

0.80 

0.40662 

0.40634 

0.00319 

0.00000 

0.90 

0.30138 

0.29967 

0.00407 

0.00000 

1.00 

0.25213 

0.24849 

0.00502 

0.00000 

1.10 

0.23991 

0.23419 

0.00597 

0.00000 

1.20 

0.25001 

0.24219 

0.00690 

0.00001 

1.30 

0.27188 

0.26197 

0.00776 

0.00001 

1.40 

0.29838 

0.28634 

0.00855 

0.00001 

1.50 

0.32489 

0.31071 

0.00925 

0.00001 

1.60  . 

0.34862 

0.33233 

0.00985 

0.00002 

1.70 

0.36803 

0.34972 

0.01036 

0.00002 

1.80 

0.38243 

0.36224 

0.01076 

0.00003 

1.90 

0.39168 

0.36982 

0.01105 

0.00003 

2.00 

0.39599 

0.37268 

0.01124 

0.00003 

2.20 

0.39156 

0.36610 

0.01130 

0.00003 

2.40 

0.373^6 

0.34675 

0.01098 

0.00003 

2.60 

0.34615 

O.3I897 

0.01035 

■ 0.00003 
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TABLE  13  continued 


Hartree-Fock  Final  wave  function 


r(a.u. 

,) 

s-type 

p-type 

d-type 

E.80 

0.31357 

0.28653 

0.00948 

0.00002 

3.00 

0.27878 

0.25237 

0.00846 

0.00002 

3.20 

0.24401 

0.21864 

0.00738 

0.00002 

3.^0 

0.21078 

0.18679 

0.00630 

0.00001 

3.60 

O.I8OO6 

0.15679 

0.00526 

0.00001 

3.80 

0.15231 

0.13177 

0.00432 

0.00001 

4.00 

0.12775 

0.10915 

0.00348 

0.00000 

4.20 

0.10635 

0.08972 

0.00277 

0.00000 

4.40 

0.08794 

0.07327 

0.00216 

0.00000 

4.60 

0.07227 

0.05949 

0.00167 

0.00000 

4.80 

O.059O8 

0.04806 

0.00127 

0.00000 

5.00 

0.04806 

0.03866 

0.00096 

0.00000 

5.50 

0.01611 

0.02207 

0.00045 

0.00000 

6.00 

0.01238 

0.00020 

0.00000 

6.50 

0.00684 

0.00008 

0.00000 

7.00 

0.00500 

0.00374 

0.00003 

0.00000 

7.50 

0.00202 

0.00001 

0.00000 

a 

C.  C 

32, 

. J.  Roothaan,  L.  M. 

186  (i960). 

Sachs  and  A.  ¥.  V/eiss,  Rev. 

Mod.  Phys 
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4,2.  PossllDle  Im-provements 

The  80  Cl  wave  function  (see  Table  10),  constitutes  a starting 
point  for  obtaining  better  approximations.  The  first  thing  we  plan 
to  do,  even  before  making  an  orbital  exponent  variation,  is  to  compute 
the  energy  of  the  80  Cl  wave  function  in  the  new  NSO  basis.  We  believe 
that  it  makes  sense  to  introduce  the  s^  and  Sg  NSO’s  from  our  final 
wave  function,  as  suggested  by  an  analysis  of  Table  14,  where  the  re- 
sults of  a c^(s^)^(s2)^  + C2(s^)^(p^)^  wave  function  are  displayed. 

In  this  table,  the  best  of  Watson's  results  shows  a relatively  high 

energy  for  the  ground  configuration,  which  is  almost  identical  with 
2 2 

the  energy  of  (s^)  (s^)  when  s^  and  s^  are  taken  to  be  the  two  NSO's 
of  highest  occupation  number  of  the  KL  shell  Cl  v/ave  function  (our 
value  is  -14.57235;  see  also  Section  3*2) . At  the  time  we  found  this 
result,  we  were  not  aware  of  Watson's,  and  very  much  in  favor  of  a 
HF  basis. 

After  this  is  done,  we  intend  to  modify  the  STO  basis  in  the 
light  of  the  results  shown  in  Table  I3.  At  the  same  time,  an  enlarge- 
ment of  the  STO  basis  is  felt  to  be  necessary.  With  a basis  of 
s,  p and  d orbitals,  it  should  be  possible  to  get  an  energy  close 

to  that  of  Kelly^^^.  Only  then  do  we  intend  to  add  higher  symmetry 

2 2 

type  orbitals.  Of  course,  configurations  like  (p)  (d)  , which  have 
been  omitted  in  this  study,  will  be  included. 

Although  this  work  is  not  completed,  it  seems  reasonable  to 
expect  that  further  significant  improvement  can  be  achieved,  and 
we  are  optimistic. 
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TABLE  14 


2 Cl  V/AVE  FUNCTIONS  FOR  Be  GROUED  STATE 


Wave  function 

(ls)^(2s)^  energy 

2 Cl  energy 

Application  of  Veillard's 
2 Cl,  SCF  formalism 

-14.572907 

-14.615449^ 

Is,  2s:  V/atson's  Hf"'^ 

: NSO  of  the  KL  shell  Cl 

-14.572986 

-14.615797^ 

Is,  2s:  HF 
2p  from  2 Cl,  SCF 

-14.57301 

-14.61594*^ 

2 Cl,  SCF 

-14.57234 

-14.61652'^ 

Clementi  and  A.  Veillard,  J.  Chem.  Phys.  44,  305I  (I966) 
E.  Watson,  Phys.  Rev.  II9,  1?0  (i960) 

^This  calculation 

R.  E.  Watson,  Ann.  Phys.  (N.Y.),  I3,  250  (I961) 
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